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ABSTRACT 

Context. Asteroseismology allows us to probe the physical conditions inside the core of red giant stars. This relies on 
the properties of the global oscillations with a mixed character that are highly sensitive to the physical properties of the 
core. However, overlapping rotational splittings and mixed-mode spacings result in complex structures in the mixed¬ 
mode pattern, which severely complicates its identification and the measurement of the asymptotic period spacing. 
Aims. This work aims at disentangling the rotational splittings from the mixed-mode spacings, in order to open the 
way to a fully automated analysis of large data sets. 

Methods. An analytical development of the mixed-mode asymptotic expansion is used to derive the period spacing 
between two consecutive mixed modes. The echelle diagrams constructed with the appropriately stretched periods are 
used to exhibit the structure of the gravity modes and of the rotational splittings. 

Results. We propose a new view on the mixed-mode oscillation pattern based on corrected periods, called stretched 
periods, that mimic the evenly spaced gravity-mode pattern. This provides a direct understanding of all oscillation 
components, even in the case of rapid rotation. The measurement of the asymptotic period spacing and the signature 
of the structural glitches on mixed modes are then made easy. 

Conclusions. This work opens the possibility to derive all seismic global parameters in an automated way, including 
the identihcation of the different rotational multiplets and the measurement of the rotational splitting, even when this 
splitting is significantly larger than the period spacing. Revealing buoyancy glitches provides a detailed view on the 
radiative core. 
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1. Introduction 


Asteroseismic observations obtained with CoRoT and 
Kepler have delivered m unificent information on the stel¬ 
lar in terior structure ('e.g.. lMichel et al.ll2008l: Caplin et ^ 
201 ill especialR for red giants ('e.g.. lDe Bidder et al.ll2009l : 


Bedding et al.H20fnll . Independent of modelling, global seis¬ 


mic properties provid e relevant estimates of the stellar 
mass es and radii (e.g., iKallinger et al.ll^OlOl: iMosser et al.l 
1201011 . The dipole modes, with a mixed character, probe 
the core and test the evolutionary sta tus of the stars 
(jBedding et al.l[2?)TTI : iMosser et al.ll2011a^ . They also pro¬ 
vide the measurement of the asymptotic period spac¬ 
ing (IMosser et al. 2012bll which i s directly related to the 
core mass ( Montalban et al]l20Ulll . Their observation also 
gives access to the dif ferential-rotation profile in red gi¬ 
ants (|Beck et al.l 1201211 . The measurement of the mean 
core ro tation for about 300 stars analyzed by_lMosser_et_al 
(l2012cll . extended toward a few subgiants (iDehenvels et al 
l2014f) . indicates that angular momentum must be effi- 


cientl y transferred from t he stellar core to the envelope 
('e.g.. iMaraues et allI2013H . At this stage, explaining the 
spinning-down r emains difficult, but a recent work by 
iBelkacem et al.l (|2015bO shows that mixed modes likely 


participate to the angular-momentum transfer and induce 
the slowing down of the core rotation. 

All these results illustrate the capability of the mixed¬ 
mode oscillations to probe the innermost radiative core and 
deliver unique information on the physical conditions deep 
inside the star. This emphasizes the need of new obser¬ 
vational constraints, on data sets as large as possible. As 
more than 13 000 red giant oscillati on spectra observe d by 
Kepler show solar-like oscillations (IStello et al.ll2013ll . ro¬ 
tation could be measured in a much larger data set than 
previously done (IMosser et al.ll201^^ . This would enable 
us to test in detail how the mean core rotation and the 
angular momentum transfer vary with stellar evolution. 

Up to now, precise rotational splittings were mostly 
manually determined. Their automated measurement is 
possible only when the rotational frequency is much smaller 
than the mixed-mode spacing. This condition is most often 
met for clump stars, but barely on the red giant branch 
(RGB). We seek a rapid, effective, and automated method 
for measuring them in thousands of stars. This first requires 
a way to disentangle rotational splittings and period spac¬ 
ings, then an automated measurement of the as ymptotic pe¬ 
riod sp acing. The method recently proposed bv iDatta et al.l 
(12015!) provides a partial answer for measuring the period 
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spacing, applicable only to non-rotating stars or stars ob¬ 
served pole-on. 


As is well known, rapid interior structure variation af¬ 
fect the sound-speed and Brunt-Vaisala frequency profiles, 
so that they modify the regular arrangement of both pure 
pressure and pure gravity modes, hence of the mixed-mode 
pattern. Such sound-spe ed glitches are clearly identified 
in red giant oscillations (|Miglio et al.l 120101: ICorsaro et al.l 
l2015t IVrard et al.ll2015all; they arise from the h elium sec¬ 
ond ionization zone fe.g.. iBroomhall et al.l[2?)l^ . Their ef¬ 
fect on the radial oscillation modes remains limited to mi¬ 
nor frequency changes compared to the large separation, 
but is enough to induce a clear difference between stars on 
the RGB or in the red clump. Thi s differ ence was observa- 
tionally sh own by Ka]lmger_et_a]Jjl20^ and theoretically 
studied by IChristen sen-Dals gaard et ^ (l2014ll . A system¬ 
atic survey by Vrard et al.l (l2015a[l firmly assesses this dif¬ 
ference and shows that it is created by acoustic glitches in 
the external envelope, related to the second ionization of 
helium. 


Buoyancy glitches are also expected in red giants, as 
a consequence of an interior structure gradient due to a 
local phenomenon related to nuclear burning or mixing 
(|Cunha et al.l 1201511 . Up t o now, they were no t observed. 
However, as suggested by ICunha et al.l (l2015ll , they may 
explain cases where the identification of the mixed-mode 
pattern is difficult or not possible. 

In this work, with a new expression of the second-order 
asymptotic expansion for dipole mixed modes, we analyze 
the structure of their oscillation pattern including rotation, 
sound-speed glitches, or buoyancy glitches. In Section[21 we 
show that the period spacing and the rotational splitting 
are constructed on similar patterns: both are bumped near 
the pressure domin ated modes. The rotatio nal pattern was 
formerly derived bv iDeheuvels et aP (l2015ll . Consequences 
of this new result are discussed in Section [31 where we show 
how the structure of the pure gravity-mode pattern can be 
revealed from the mixed-mode pattern. Section|3]shows how 
rotation can be fully disentangled from the mixed mode 
pattern, even in the most complex cases where rotational 
splittings and period s pacings largely overl ap (see for in¬ 
stance Figs. A.4-A.6 of iMosser et al.lF2012cll . We then in¬ 
vestigate the principle of the full automation of the analysis 
of the red giant oscillation spectrum and test the influence 
of structure glitches in Section [31 

A companion paper, IVrard et~aI1 (|2015bll . presents the 
automated method for deriving the asymptotic period spac¬ 
ing in Kepler red giants, which is based on the present work. 


2. Mixed mode period spacing 


Follow ing the work of IShibahashil (Il979ll and lUnno et al.l 
(|1989ll . asymptotic expansions of mixed modes have 
been derived for different seismic parameters: eigen- 
frequencies (iMosser et al 201291 . period spacings 
( Christensen-Dalsgaardl 2012ll. rotati onal splittings 


( Goupil et al.ll201^ Deheuvels et al.ll2015ll . 


Here, we intend to provide expansions taking into ac¬ 
count complex expressions of the pure pressure and pure 
gravity contributions. 


2.1. Asymptotic expansion 

IShibahasliil (|l979ll and lUnno etaTI (Il989ll derived an im¬ 
plicit asymptotic relation for mixed modes, which expresses 
as 


tan^p — gtan^g, 


(1) 


where q is the coupling factor between the gravity and pres¬ 
sure components of the modes and where the phase Op and 
Og refer, respectively, to the pressure- and gravity-wave con¬ 
tributions. At first order, the phases for dipole modes are 
related to the period spacing AHi and large separation Av, 
respectively. They write 

2 ’ 

as used for red giants by IMosser et al.l (|2012bl) . but here 
correctly accounting for the contribution £/2, with i = 1, 
and for the offsets Eg and £p. 

Owing to the form of the mixed-mode asymptotic ex¬ 
pansion, we can modify the phases in Eqs. m and ([3]) with 
additional terms, provided their contributions are multi¬ 
ples of TT. This is a useful artifice for introducing the eigen- 
frequencies of either pure pressure or pure gravity modes. 
Therefore, we rewrite the phases as 


( 2 ) 

( 3 ) 


TT-^ f- - — V 

AHi \v Vg) 

V -Vp 

Av{npy 


( 4 ) 

( 5 ) 


where Vp and Vg are the asymptotic frequencies of pure 
pressure and gravity modes, respectively, and Az/(np) is the 
frequency difference between the consecutive pure pressure 
radial modes with radial orders rip and rip -|- 1. The offsets 
and 112 term introduced in Eqs. (I2I3I1 are included in the 
terms Vp and Vg. This means that we now have the pos¬ 
sibility to use Eqs. (gH) at any order of the asymptotic 
expansions for the pure pressure and gravity contributions. 

For iTg, acceptable fits of red giant oscillation spectra 
are b ased on the first-order asymptotic expansion (iTassoull 
1198011 . For dipole modes, we have 


^ = (-ng-keg)Ani, 

i/g 


( 6 ) 


where rig is the gravity radial order, usually defined with a 
negative value, and the offset Eg = 1/4 — Eg'. The gravity 
offset Eg' is a small but complicate function sensitive to 
the stratification near the boundar y between the radiative 
core and the convective envelope (iProvost fc BerthomieuI 
Il986ll . If desirable, we may consider in Eq. JIL a second- 
order expansion for the gravity expansion (ITassoull 1198011 . 
or the contribution of bu oyancy glitches fe.g.. iMiglio et al.l 
120081 : [ciinha et ai]l2015ll . This is done in Section |5l 

For Vp, the high quality of the seismic dat a requires the 
use o f the second-order asymptotic expansion (IMosser et al.l 
1201311 . For dipole modes, we have 


^(^p) — 


1 a . 

2 + + "2 ~ 


y + doi ) Av, (7) 
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Fig. 1. Relative period spacings AP/AIIi and relative ro¬ 
tational splittings C for a synthetic RGB mixed-mode spec¬ 
trum. The red line shows the asymptotic function; the 
blue crosses represent the period spacings derived from the 
asymptotic expansion of mixed modes and scaled to AIIi. 
The dashed line shows the minimum values Cmin reached by 
pressure-dominated mixed modes; the dotted line shows the 
maximum values Cmax reached by gravity-dominated mixed 
modes. 

where Up is the pressure radial order, a represents the cur¬ 
vature of the radial oscillation pattern, rimax = ^'max/Az/ — 
Ep is the non-integer order at the frequency z^max of max¬ 
imum oscillation signal, and doi is the small separation, 
namely the the distance, in units of Az.^, of the pure pres¬ 
sure dipole mode compar ed to the midpoint between the 
surrounding radial modes (iMosser et al.ll201^ . In this case, 
the observed large separation Az/(np) in Eq. ([5]) increases 
with increasing radial order: 

Az/(np) = Az/ (1-I-a(np - Umax)), (8) 

which is obtained from the derivation of Eq. ©. As for 
the gravity contribution, the asymptotic expressions of the 
pure p modes can be as precise as required; Vp may include 
the g litch component due to the second ionization of he- 


hum 1 

Midio et alJl2010HBroomhall et al.ll20l4lVrard et alJ 

I2m5a 

). 


2.2. Bumped period spacing 


Fig. 2. Comparison of the function C considered as a con¬ 
tinuous function of frequency (continuous line with a color 
modulated by the value of cos^ 0g) with the values de¬ 
rived for the mixed-mode frequencies (full squares with a 
color modulated by the value of cos^0p). Values of cos^ 
close to 0 (1) appear in light blue (dark blue). For gravity- 
dominated mixed modes, the phase term 9p acts for having 
(C values close to the local maxima; for pressure-dominated 
mixed modes, the C values are close to the local minima 
Cmin (dashed line). The simulation is based on the second- 
order asy mptotic press u re-mo de pattern proposed for red 
giants bv iMosser et al.l (l2013l) . with Ai/ = 5/rHz, and on 
the hrst-order asymptotic expansion of gravity modes, with 
AIIi = 65 s and Eg = 0. 


when the variation of Az/(np) with frequency is neglected 
in the pressure term. One easily obtains 


2 dP Av(np) cos^ Op 

V - = q - — -^ 

drim AIIi cos^ 0g 


dP 

dUm 



hence 

1 dP r. I i^^AIIi cos^0g 

-= 1-1 -^ 

AIIi dum [ q Av{np) cos^ Op 


( 11 ) 


( 12 ) 


The case drim = 1 corresponds to the period difference 
between two consecutive mixed modes, and thus provides 
us with the relative bumped period spacing 


In order to retrieve the relative period spacing AP/AIIi 
between two consecutive mixed modes, one needs to intro¬ 
duce the mixed-mode radial order, = rip + rig. We then 
choose to write the period P = 1/z/ of a mixed mode as 


AP 

ABi 




where the definition 


(13) 


P = UmAIIi -I- p, (9) 

where the residual term p expresses the departure to an 
evenly-spaced comb function. This period is introduced in 
the phases Og and Op provided by Eqs. m and and 
linked by the asymptotic relation (Eq. [T]). In the gravity 
contribution, all terms multiple of tt cancel, owing to the 
property of the tangent function, so that with Eq. m the 
term tan0g reduces to tan 7r(p/AIIi ~ Sg)- The derivation 
of Eq. ([T|) with respect to rim then gives 

d / 1 \ dp 

drim yrZmAIIi +pj _ drim 

Az^(np) cos^ 0p ^Allicos^ffg’ 


C(^) = 


1 z^^ABi 


COS^ TT 


ABi 


q Av{np) 


Au{np) 


(14) 


matches the first-orde r asymptotic expansion derived by 
iDeheuvels et al.l (|2015f) for expressing the mixed-mode ro¬ 
tational splitting 


^^rot — ^^rot,g C ^U^ot p (1 ^), (1^) 

where Surot.g and Surot^p are the rotational splittings related 
to pure gravity or pure pressure modes. In the limit case 
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where the mean envelope rotation is negligible compared to 
the mean core rotation, this simplifies into 




^^rot ,g 


(16) 


Hence, we have demonstrated that the period spacing com¬ 
pared to AHi (Eq. [TS]) and the rotational splitting com¬ 
pared to (5r'rot,g lEg. 1161) follow the same distribution. This 
result emphasizes that the coupling between the pressure 
and gravity terms plays the same role for arranging the 
mixed-mode peri ods and the ( f reque ncy) rotational split¬ 
tings. As shown in lGoupil et al.l (l2013ll and iDeheuvels et ^ 
(|201fill . this behavior is governed by the inertia of the mixed 
modes: Q is the ratio between the kinetic energy in the 
radiative cavity and the total kinetic energy. The equal¬ 
ity between AP/AHi and C is seen in Fig. [I] where the 
bumped period spacings AP defined as P{n^ + 1) ~ P{nm) 
are compared to the function (^. In order to avoid quanti¬ 
zation biases, we have plotted AP/AHi at the abscissae 
(P(nn,) + P(nn, + l))/2. 

We must stress that the values of C, obtained for the dis¬ 
crete mixed-mode frequencies v do not reflect the complex¬ 
ity of the function (( if considered as a continuous function 
of frequency. For clarity, we denote the discrete values 
and Cif) the continuous function. The contributions of the 
phases 9p and 9g to C(/) are shown in Fig. O 

Interestingly, the global pattern (^{v) reached for mixed 
modes mainly depends on the asymptotic parameters of 
the pure pressure modes, since the location of the local 
minima only depends on the frequency of the pure pres¬ 
sure d ipole modes. Conve rsely, it hardly depends on AHi 
and q (IVrard et al.ll20I5bll . Limiting cases corresponding to 
gravity-dominated mixed modes (g-m mode) and pressure- 
dominated mixed modes (p-m mode) are derived in the 
Appendix. 





Fig. 3. Interpolation of the function ( for the period spac¬ 
ing AHi, derived from the values obtained with various 
periods in the range AHi)! ± r'i„axAni/2), where t'max is 
the frequency of maximum oscillation signal. 


By comparison with Eq. ([T}, we get $ = 7r/2 — 0g. The 
term d$/di/ derived from Eq. (1191) ensures the agreement 
between the different expressions of AP/AHi. We note that 
the use of Eq. (HI avoids the introduction of the derivative 
term present in Eq. m- As rapid variations arise from the 
9g contribution, our new expression of AP/AHi (Eqs. [131 
lUl) provides a more precise numerical result than Eq. dm 
requiring numerical derivative. 

3. Stretched periods 

In this section, we show how the function (( can be derived 
from the analysis of the pressure mode pattern and how 
it can be used to extract information on the pure gravity 
mode pattern. 


2.3. Comparison with previous work 

Recently, iBenomar et al.l (l2014f l have remarked that the ra¬ 
tio Ii/Iq of the inertia of dipole modes compared to radial 
modes may provide an alternative to using mixed-mode fre¬ 
quencies and give essentially the same information as the 
mixed-mode asymptotic equation. We confirm that this is 
the case: their function Ii/Iq can be related to ( if one as¬ 
sumes that the difference between the total inertia of mixed 
modes and the contribution of the radiative cavity only 
exactly corresponds to the inertia of radial modes. Their 
equation (10), with the factor <7^/4 corrected into q, can be 
written 


h 

h 


1-C 


(17) 


A different expression of the bumped period s pacin g 
was earlier introduced by IChristensen-DalsgaardI (l2012ll . 
Rewriting Eq. (34) of his work with our notation gives 


AP 

AHi 



z/^AHi d$' 
TT dlV 


(18) 


with $ defined in IChristensen-DalsgaardI (|2012f) by 


tan $ =-—. 

tan 9p 


(19) 


3.1. Identification of the mixed modes 

In a red giant oscillation spectrum, the location of 
the mixed modes is in fact fixed by the measurement 
of the large separation Av. The determination of Az/, 
first derived from the envelo pe autocorrelation function 
(jMosser fc Appourchauxl [200911 . then refined with the un i- 
versal red giant oscillation pattern (iMosser et al.ll20lTbll . 
provides the identification of the radial modes and helps 
locating the frequency ranges where mixed modes cannot 
be mistaken for radial or quad rupole modes . The b ack- 
ground parameters, derived as in iMosser et al.l (l2012all . are 
used for correcting the granulation contribution in the fre¬ 
quency range around Umax where oscillations are observed. 
Hence, mixed modes can be automatically identified in fre¬ 
quency ranges having no radial and quadrupole modes, with 
a height significantly above the background. 

3.2. Estimate of the function ( 

The same information is used for deriving the function 
The precise location of the minima of the function is ob¬ 
tained from the location of the p-m modes, which is fixed by 
the radial pa ttern. Then, basic in formation from the AHi - 
Au diagram (jMosser et al.l[20T^ is enough to obtain a pre¬ 
cise estimate of (. Large modifications of AHi only slightly 
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Fig. 4. Stretched period r as a function of frequency (dark 
blue line), compared to a reference varying asl/v (dashed 
line). The derivative dr/dv is superimposed (light blue 
line), with an arbitrary scaling factor to fit in the same 
window as r; the dotted line is representative of the term 

l/i/2. 

modify the C profile (IVrard et al.ll20159l . Furthermore, the 
variation of q with AIIi, hence with Av, also obeys a deter¬ 
ministic relation, so that q can be predicted with a precision 
better than 20 %, which is enough. 

As shown in Fig. [2j the function is properly defined 
for the mixed-mode frequencies only. Therefore, to obtain 
a continuous function representative of the period change 
due to the mixing of the modes, we have to interpolate the 
values of C(*^)- Figure |3] graphically explains how this can 
be obtained with a small modulation of the period spac¬ 
ing, assuming that ( only sli ghtly depends on A IIi . This 
property is fully developed in lVrard et al.l (l2015bll . 

3.3. Stretched periods 

The interpolated form of C is then used to turn the frequen¬ 
cies into periods r with the differential equation 


from which we can integrate the periods t of mixed modes 
(Fig. 2]). The constant of integration can be arbitrarily fixed 
since it plays no role. Owing to the ( profile (Fig.[T|), we call 
them stretched periods. They are used for drawing echelle 
diagrams, where v is plotted as a function of r modulo AIIi. 


period used for folding the stretched period influences the 
slope of the ridge. It is easy to ensure a vertical alignment 
and then extract the measurement of the pe riod spacing 
AIIi. This is done in the companion paper (I Vrard et al.l 
l2015bD and we let the discussion on the performance for 
measuring AIIi to this work. 

3.5. Precise measurement of the coupling factor 

The precision on the measurement of the coupling term q 
can be significantly enhanced. Previous estimates showed 
unc ertainties as large as 0.05 for a mean value of q about 
0.2 (|Mosser et al.ll201^ . Here, the expression of Cmin gives 
a new way for measuring q in an automated manner. As 
in Figure [TJ one can compare the period spacings, for a 
fixed m ridge, with the function CminAHi, and derive q 
with a smaller uncertainty. It is even possible to reduce 
this uncertainty to about 0.02 by tuning q for ensuring a 
minimum spread around p-m modes in the echelle diagram 
constructed with the stretched periods. 

4. Disentangling rotation 

In this section, we intend to use the analytical expression 
of the period spacing (Eq. [T3]) to disentangle the complex 
structure of the mixed mode pattern, especially for red gi- 
ants on the RGB wh ere the most complex cases are met 
(iMosser et al.l l2012cD . For this purpose , we have consid¬ 
ered st ars analyzed in iBeck et al.l (|2ni2ll and iMosser et al.l 
(|2 oT^I3i . We have retrieved the public data set in the 
Kepler archive in order to benefit from the longest avail¬ 
able time series (about 4 years). 

4.1. Rotational splittings 

The rotation being considered as a perturbation of the 
dipole mixed-mode pattern and the surface rotation be¬ 
ing negligible in most cases, the unpe rturbed frequency v 
is changed by rotation into u — mdurot (|CIouDil et al.ll201.^ . 
According to Eq. (fT5|) and to the profile of C, this reduces 
to V — mC(5;zrot,core, cven in the cases where the envelope 
rotation canno t be neglected, as for instance in the sec¬ 
ondary clump (iDeheuvels et al.l[2015ll . For dipole modes, 
azimuthal orders m S {—1,0,1}. Hence, the differences in 
the stretched periods of dipole modes with same azimuthal 
order m can be approximated by 

Ar^^AHi + (21) 


3.4. Echelle diagrams based on the stretched periods 

We first consider the case of a star seen pole-on, for 
which rotation is not an issue; the full case including ro¬ 
tation is solved in Section m The left panel of Fig. [S] 
shows the varying period spacings of the mixed modes. 
The echelle diagram based on the periods l/u (middle 
panel) shows the c lassical S-shape (iBedding et al.l 120111 : 
IMosser et al.ll2012bll . whereas the echelle spectrum based 
on the stretched periods r (right panel) exhibits a single 
nearly-vertical ridge, where discrepant peaks can be easily 
identified as speckle structures due to the short-lived p-m 
modes or as £ = 3 modes. As for classical echelle diagrams 
used for estimating the large separation Az/, the value of the 


if w e assume that the m ean envelope rotation can be omit¬ 
ted (|GouDil et al.ll201^ . Even for an important rotational 
splitting, as large as or even larger than the period spacing, 
the relative correction due to rotation is small. Therefore, 
the stretched period spacings Arm=±i are close to AHi 
(Table [TJ . Since modes are observed in a limited frequency 
range around the frequency zzmax of maximum oscillation 
signal, we can approximate the varying term (^/u of Eq. (1211) 
by N/{{N -I- l)t'max), where the function 


N{u) 


Av 

Aniz^2 


( 22 ) 


shows how the period spacings compare to the large sepa¬ 
ration Av. It corresponds to the number of gravity modes 
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Fig. 5. Echelle diagrams of the RGB star KIC 9882316 that shows no rotational signature. Left: Classical frequency 
echelle diagram as a function of the dimensionless frequency v/lS.v — (n-p + e). The radial order Up is indicated on the 
y-axis. Radial modes (highlighted in red) are centered on 0, quadrupole modes (highlighted in green), near —0.12 (with a 
radial order Up — 1), and £ = 3 modes, sometimes observed, (highlighted in light blue) near 0.22. Dipole mixed modes are 
identified with the frequency given by the asymptotic relation of mixed modes, in /iHz. The fit is based on peaks showing 
a height larger than eight times the mean background value (grey dashed lines). Middle: Classical period echelle diagram, 
where the abscissa is the period modulo the period spacing AIIi. The most prominent mixed modes, marked with blue 
filled squares (in light blue for peaks in the vicinity of p-m modes), are automatically identified. In the background of the 
figure, the spectra are plotted twice and top to tail for making the mode identification easier, with the pressure radial 
orders indicated on the radial modes. Right: Stretched period echelle diagram, where the abscissa is the stretched period 
T modulo the period spacing AIIi. Pressure-dominated mixed modes are coded in light blue. 


Table 1. Rotational splittings 


Pressure modes Hand-made analysis Analysis based on stretched periods 


KIG 

i^max 

/iHz 

Av 

/iHz 

AHi 

s 

A/* (/^max) 

rot,core 

iiHz 

ATtti—1 

s 

ATjti —_1 

S 

(Ar^n) 

s 

^rot 

% 

rot,core 

nHz 

6144777 

128.4 

11.029 

79.0 

8.5 

242 

78.75 

79.25 

79.0 

0.34 

227 

9267654 

118.0 

10.358 

78.0 

9.5 

490 

77.30 

78.55 

77.9 

0.75 

523 

12008916 

160.8 

12.906 

79.8 

6.3 

430 

79.45 

80.15 

79.8 

0.46 

409 


fal rotational splitting derived from the use of lGoupil et al.l (1201,11 1. 
(b) rotational splitting derived from this work 


per Ai^-wide frequency range around the frequency iz. This 
shows that the period spacing of each component of the 
rotation multiplet is nearly uniform, close to 

Arm ~ ABi (1 -I- mxrot), (23) 

with 

A/'fz^max) ^^rot.core 

2 777-TTT-■ 

Af (i^max) T 1 ^max 

We note that the ratio Af/{M -b 1) is a consequence of the 
addition of an extra pressure mode among the Af gravity 
modes that construct the {Af + 1) mixed modes in a Az/- 
wide frequency range. 


4.2. Measurement of ABi and dz/rot.core 

The use of Eq. (1^ allows us to correct the variation of 
the rotational splitting with frequency and their asymme¬ 
try near the p-m modes. According to Eq. (1^7)) . each az¬ 
imuthal order forms a well-identified ridge in the echelle 
diagram. Each ridge approximately shows linear variation 
of the rotationally perturbed period spacings (Fig.lHl), from 
which we can derive the period spacing AIIi and the rota¬ 
tional splitting dz/pot- When the rotation splitting exceeds 
half the mixed-mode spacing, the m = ±1 ridges intersect 
(Eig. [6l right panel). Results obtained for these stars are 
given in Table [TJ Values agree with th e rotational split¬ 
tings measured in iMosser et al.l (l2012cll , with a small dif¬ 
ference due to the approximation used in this earlier paper, 
which is explained by the difference between the approxi- 
mate form of the rotational splitting found by (jCoupil et al.l 
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KIC1200B916 - 79.80 s 



T mod An^ (s) 


KIC 6144777 - 79.20 s 



T mod Aril (^) 


KIC 9267654 - ABj^ 77.90 s 



T mod Aril (^) 


Fig. 6. Echelle diagrams based on the stretched periods of dipole mixed modes. The abscissae are the periods r, derived 
from Eq. (uni), modulo the asymptotic period spacing AIIi; the ordinates are frequencies. The most prominent mixed 
modes, marked with filled squares, are automatically identified. In the background of the hgure, the spectra are plotted 
twice and top to tail for making the mode identification easier, with the pressure radial orders indicated on the radial 
modes. From the left to the right, we show RGB stars with different rotational splittings. Left: Star seen equator-on, 
where only m = ±1 modes are visible, and with a moderate mean core rotation; they draw two ridges from which the 
rotation frequency Sv^ot can be derived. Middle: Star seen with an inclination such that all three components of the 
rotational multiplets are visible, and with a moderate mean core rotation. Right: Star seen equator-on; as a consequence 
of the large mean core rotation, the ridges m = ±1 intersect. 


l2013ll and the more pr ecise asymptotic rotational splitting 
( Deheuvels et al.]l2015ll . 

When the surface rot ation cannot be omit t ed, a s is 
the case on the low RG (iDeheuve ls et al.l 201^ 2 01411 or 
in the secondary red clump ( Deheuvels et al.l 2015[l . the 
rotation correction C^^^rot in Eq. (|2T) must be change d 
into C(<J^'rot ,core Sl^rot ,env ) ^^rot ,env iGoupil et al.l[ 201 ^ . 
Gonclusions remain the same, but estimating the separated 
contribution of the core and of the envelope then requires 
a separated study of p-m and g-m modes. 


4.3. Toward an automated determination of the global 
seismic parameters 

The method presented above has shown its ability to ex¬ 
hibit in a clear way the rotational structure of the mixed 
modes and identify their azimuthal orders. This comes from 
the fact that Eq. (EOl) allows us to directly analyze the 
period pattern, and helps us avoiding the identification of 
the mixed modes pri or to the measure ment of AIIi as in 
other methods (e.g., iDatta et al.ll2ni.^ . In fact, in many 
cases and especially on the RGB, defining period spacings 
of dipole mixed modes requires first the full identification of 
the dipole modes, which is a severe drawback for complex 
spectra: rotational splittings, £ = 3 modes, and short-lived 
p-m mode conspire against the identification of the angular 


degree and azimuthal order of the modes, then against the 
measurement of relevant AP. The use of Eq. (uni) alleviates 
this problem. 

All stages used to draw the echelle diagrams were based 
on automated methods. The steps for obtaining Av were 
recalled in Section [3] and are not modified by rotation. The 
asymptotic peri od spacings can be obtained in an auto¬ 
mated way too (|Vrard et al.l l2ni5'9l . The next step is to 
make profit of the views developed in this work to automate 
the measurement of the rotational splitting too; this is the 
aim of a forthcoming paper. The met hod can also be used t o 
investigate possible rapid rotators (lOuazzani et al. I [201,311 . 
Possible complications arising from structural glitches are 
presented in the next Section. 

5. Glitches 

In this Section, we analyze how structural glitches modify 
the mixed-mode pattern and might complicate the oscilla¬ 
tion spectrum. Glitches are due to rapid variation either in 
the sound-speed profile or in the Brunt-Vaisala profile. 

5.1. Sound-speed glitches 

The pressure mode spectrum of red giants is modulated by 
acoustic glitches due to abrupt variation of the first adi- 
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abatic exp onent at the locat i on of the second ioniz ation 
of helium (iMiglio et alJ l201fll : iBroomhall et al l l2014li. We 
have m odelled the pressure glitches following IVrard et alJ 
(|2015ari . with a simple sinusoidal modulation added to the 
pure pressure modes Vp. We considered two cases, corre¬ 
sponding to typical clump and RGB stars, and used the 
mean values of the glitch periods and amplitudes given in 
IVrard et alJ (|2015aD : the periods of the glitch signature are 
in the range [3-4], in units Az/, and the amplitude about 
a few percents, in units Az/ too. We note that the introduc¬ 
tion of the sound-speed glitch does not modify the scheme 
presented in Section [21 so that one can easily retrieve the 
period spacings (Fig. jT]). 

The analysis of the stretched periods confirms that only 
p-m modes are significantly affected by the sound-speed 
glitches. Their location, hence the value doi of the small 
separation in Eq. must be precisely determined in order 
to perform a correct analysis of the glitch signature. 

The influence on the g-m modes is reduced since the 
pure gravity-mode pattern is not affected by the sound- 
speed glitch. According to the amplitude observed for 
acoustic glitches, the deviation in the p-m mode is limited 
to a small fraction of AIIi. The relative amplitude Xp of 
the acoustic glitch, which corresponds to a frequency shift 
of XpAiy, translates into a period shift of the pure pressure 
modes about Xp M AIIi. The maximum period shift Sxp of 
the mixed modes is much smaller, since the gravity compo¬ 
nent of the wave is not changed. We provide an estimate 
of this effect with a set of simulation of synthetic glitches, 
performed at various evolutionary stages. Scaling relations 
of 5tp with Az/, z.^max and AIIi can be summarized by 

STp^OAxp A/'(z/„,ax) ABi . (25) 

Obser ved values of Xp/Av are about 1 - 2 % (|Vrard et al.1 
l20I5all . and at there are about 10 gravity modes per 
Av frequency range, so that the perturbation ^Tp/AIIi is 
at most about 1 - 2%. In any case, sound-speed glitches 
are neither expected to perturb the determination of the 
mixed-mode pattern, nor the measurement of AIIi ■ 


5.2. Buoyancy glitches 

For bu oyancy glitches, we use the analysis of iMiglio et al.l 
(|200i^ describing them in SPB and 7 Doradus stars, since 
the analytical expressions they derive are general enough 
to be applied to the red giant case. The asymptotic period 
spacing, modulated by a buoyancy glitch, can be written 

Anigiii,h = ABi (l + 5 (l/zz)), (26) 

where g is a periodic function. Its period is defined by 


g 


ABi 





(27) 


where Nbw is the Brunt-Vaisala frequency, r^i {ri\[ 2 ) is the 
inner (outer) boundary of the radiative core, and where Cg 
is the radius of the buoyancy glitclQ. 

^ We note that this definition implies that the ratio Q is neces¬ 
sarily larger than unity. It is also possible to change the bound¬ 
aries of the denominator into, respectively, rjvi and rg, for an¬ 
other definition Q* of the period, which verifies 1/C/* = 1 — 1/C/. 




Fig. 7. Stretched period echelle diagrams of synthetic 
mixed-mode spectra including sound-speed glitches. Left: 
RBG star; Right: red clump star. Horizontal dotted lines 
indicate the location of the pure pressure dipole modes. 



Fig. 8 . Relative period spacing AP/AHi (blue line) de¬ 
rived from the function ( calculated without glitches (red 
line and plusses). Both values of C and C — dg/dn^ are 
superimposed and cannot be distinguished. 


When Brunt-Vaisala glitches are present, the derivation 
of Eq. (m for obtaining Eq. (fldll must take the variation of 
g into account. We suppose that the buoyancy glitch can 
be treated as a perturbation, so that the influence of the 
glitch in the phase 9p can be neglected. Hence, Eq. m is 
changed into 


1 dP dg 

-= C- — 

AHi dum drim 


(28) 


The correcting term on the right-side of this equation in¬ 
cludes the glitch contribution. We have checked that, except 
for huge glitches, the glitch-perturbed correction (/g is very 
close to the glitch-free function ( (Fig.jS]). As a consequence, 
the function ^ can be safely used in Eq. (EHl) instead of Q to 
derive stretched periods. The stretched periods calculated 
with C exhibit the glitch signature; conversely, when Cg is 
used in Eq. (EHl), the glitch signature is corrected. We also 
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Ani = 80 s; A = 0.95; G=1.01 Ani = 80 s; A = 0.95; G=l.l Ani = 80 s; A=0.95; G = 2.0 Ani = 80 s; A = 0.95; G=11.0 Ani = 80 s; A = 0.95; G=101 
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Fig. 9. Stretched period echelle diagrams of synthetic mixed-mode spectra including buoyancy glitches. The period of 
the buoyancy-glitch modulation is indicated for each plot: it increases from the left to the right. All other parameters 
are the same in all plots. Horizontal dotted lines indicate the location of the pure pressure dipole modes. 



notice that the expression used for the glitch is coherent 
with the modelling done bv ICunha et al.ll2015ll . 

5.3. Deep and shallow buoyancy glitches 

We analyzed the cases where the function g expressing the 
buoyancy glitch varies sinusoidally and tested glitch periods 
Q in a broad range, from 1.01 to 100. Even if the sine form 
is a simplification, it helps investigating the glitch signature 
since it carries the most important information, namely the 
location of the glitch, which is related to the re lative period 
Q. We denote the amplitude of the glitch as in iMiglio et al.l 
( 200811 : the amplitude A corresponds to a relative jump 
(1 — A) in the Brunt-Vaisala profile at the glitch location. 

Figure ini shows that the echelle diagrams based on the 
stretched period spacings make it possible to retrieve the 
glitch profile. We note that the glitches with relative peri¬ 
ods {Q — 1) and l/(^ — 1) show similar modulation, except 
near the p-m modes. Such a degeneracy in the glitch in¬ 
formation is well known: for symmetry reasons, deep and 
shallow glitches cannot be distinguished. Here, the degen¬ 
eracy is broken by the p-m modes. We might distinguish 
three main cases: 

- For Q values close to unity, that is for glitches located 
near the lower boundary of the Brunt-Vaisala cavity, the 
period bumping near the p-m modes induces a significant 
shift compared to the period ^AHi (Eq. [2S1). As a result, 
the large variation of the phase of the modulation exceeds 
27r so that it cannot be corrected. So, a rapid variation of 
the stretched period near the p-m modes is the signature of 
a shallow buoyancy glitch. The amplitude of the spikes at 
the p-m modes can be used to investigate the amplitude of 
the modulation of the glitch. According to Fig.O the spikes 
near the p-m modes have the same amplitude as the global 
modulation. This result is robust; it derives from a phase 
effect and is independent of the exact form of the periodic 
function g. 

- For periods Q larger than 2, the shift at the p-m mode 
remains limited, so that the C, contribution makes it pos¬ 
sible to retrieve the glitch function g{l/v) in the echelle 
diagram based on the stretched period. Large values of Q, 
corresponding to shallow glitches, show long-period modu¬ 
lations. The absence of any accident near the p-m mode in 


KIC 9332840 - Ani = 301.00 s 



-150 -100 -50 0 50 100 150 

T mod AIlj (s) 


Fig. 10. Stretched period echelle diagram of KIC 9332840, 
showing a large-period modulation (dashed red line) com¬ 
patible with a shallow buoyancy glitch. Peaks are identified 
as in Figs. [S] and [HI The minimum and maximum values of 
the varying period are indicated in orange. 


the echelle diagram is the signature of a shallow glitch. 

- The case 1/ = 2 is an intermediate case: the glitch sig¬ 
nature looks like a sawtooth profile; the varying amplitude 
of the modulation is a moire effect. Determining the mean 
value of AHi may then be difficult. However, such cases 
were not encounter ed in the large data set analyzed in 
iMosser et al.l (l2014ll . 
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5.4. An example of buoyancy glitch 

We consider the case of the red giant star KIC 9332840 
observed by Kepler. This star has reached t he red clump 
and sh ows a complex mixed-mode pattern (iMosser et al.l 
l2012all . Two values of the period spacing were reported, de¬ 
pending on the frequency range: 298 s in the low-frequency 
range and 306 s in the high-frequency range. This variation 
can be attributed to the modulation induced by a buoy¬ 
ancy glitch (Fig. [ini). An alternative explanation should 
be a deviation to the first-order asymptotic expansion. 
Observationally, the observed mixed modes have so high 
gravity radial orders, about 80, that we do not favor this 
hypothesis. Theoretical simulations bv ICunha et al.l (|2015D 
seem to confirm this for RGB stars. 

When the modulation is interpreted as a buoyancy 
glitch, the stretched period diagram of this spectrum shows 
a signature with a high Q value, in the range [40 - 50]. This 
means that the glitch has a shallow location and might be 
related to some accident in the Brunt-Vaisala profile near 
the outer boundary of the dense radiative core. The ampli¬ 
tude A, about 0.08±0.02, provides a measure of the glitch 
contrast in the Brunt-Vaisala profile. 

So, contrary to sound-speed glitches, buoyancy glitches 
may significantly perturb the mixed-m ode pattern. 
Howe ver, according to previous observations (iMosser et al.l 
l2014fl . this only occurs in a small minority of cases. A com¬ 
prehensive observational study of buoyancy glitches further 
requires a systematic analysis of the red giant oscillation 
spectrum. 

6. Conclusion 

We have rewritten the asymptotic expansion of dipole 
mixed modes in order to account for the most precise de¬ 
scription of the pure pressure and pure gravity contribu¬ 
tions. We have shown that the bumped period spacing and 
frequency rotational splitting follow the same patter n con¬ 
structed by the function ( (iDeheuvels et al.M2015ll . This 
equality implies that the information present in the rotation 
splittings and in the period spacings is degenerated: both 
signatures are derived from the ratio of the kinetic energy 
in the radiative core compared to the total kinetic energy 
of the modes. Then, we have shown how the function can 
be used for stretching the mixed-mode pattern and deriving 
the contribution of the pure gravity-mode pattern. In fact, 
stretched periods mimic the pure gravity period spacing: 
they are evenly spaced, and the mean spacing corresponds 
to the asymptotic period spacing AHi. As a result, each 
component of the rotational multiplets can be identified 
and the asymptotic period spacing can be measured even 
in presence of important rotational splittings. 

It follows that all steps of the red giant spectrum anal¬ 
ysis can be automated, since the complex mixed-mode for¬ 
est is now as regular as an artificially planted thicket. We 
have derived an a utomated method fo r measuring AHi in a 
companion paper (|Vrard et al.ll2015bll ; the automated mea¬ 
surement of the core rotation in 13 000 red giant spectra 
observed by Kepler and showing solar-like oscillations is in 
progress. 

Another output is the possible interpretation of the 
modulation of the period spacing as buoyancy glitches. 
Echelle diagrams based on the stretched periods are able 
to put them in evidence. This reinforces the capability 


of mixed modes for probing the stellar cores and infer¬ 
ring unique information on the physical conditions in the 
nuclear-burning region. 
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Appendix A: Limiting cases 

A.l. Gravity-dominated mixed modes 

The frequency v oi a, gravity-dominated mixed mode (g-m 
mode) is significantly offset compared to the pure pressure 
frequency Vp, so that its phase Op is close to tt/2, mod¬ 
ulo TT and its cosine is small. Hence, the function shows 
large variation. For these g-m modes, tan0p is large, so that 
tan0g is necessarily large too, according to the asymptotic 
relation (Eq. [T]). Considering that both cosine values are 
small, one can approximate them by the inverse of their 
tangents. As a result, when taking the asymptotic expan¬ 
sion into account, the value of ( for g-m modes is close to 
Cmax defined by 


r. 1/2 ahi 1 

-1 

r. Q 1 

1 + 9 A / 

[ Ai/(np)J 

— 



(A.l) 


hence close to unity in most cases except on the lower part 
of the red giant branch. As a by-product of this discussion, 
we note that the g-m periods are close to (rim 1/2 -\- 
eg)Ani. So, even if g-m modes behave as gravity modes, 
their frequencies are shifted due to an extra phase added 
by the coupling; this situation resembles the tt phase added 
to the wavefront of an optical ray crossing a focal point. 


A.2. Pressure-dominated mixed modes 


On the contrary, for pressure-dominated (p-m) mixed 
modes, the phase Op is close to 0: the function ( shows 
reduced variation between Cmin and 1. The low value of 
tan dp implies the low value of tandg; cos^ dp and cos^ dg 
are both close to unity, so that for p-m modes one finds ^ 
close to (Cmin defined by 


r 11/2 AHi 1 

-1 

r, 11 

1 -f . , . 

q An (Tip) 

— 



(A.2) 


This term Cmin is significantly different from unity; as is well 
known, period spacings of p-m modes can not be used to 
derive the asymptotic period spacing (e.g., iBedding et al.l 
[Ml iStello etal][20Tl . 
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